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Abstract. We report the first interferomteric detection of 1 83 GHz water emission in the low- 
mass protostar Serpens SMM1 using the Submillimeter Array with a resolution of 3" and rms 
of ~7 Jy in a 3 km s _1 bin. Due to the small size and high brightnessof more than 240 Jy/beam, 
it appears to be maser emission. In total three maser spots were detected out to ~ 700 AU from 
the central protostar, lying along the red-shifted outflow axis, outside the circumstellar disk 
but within the envelope region as evidenced by the continuum measurements. Two of the 
maser spots appear to be blue-shifted by about 1 to 2 km s _I . No extended or compact thermal 
emission from a passively heated protostellar envelope was detected with a limit of 7 Jy ( 1 6 K), 
in agreement with recent modelling efforts. We propose that the maser spots originate within 
the cavity walls due to the interaction of the outflow jet with the surrounding protostellar 
envelope. Hydrodynamical models predict that such regions can be dense and warm enough 
to invert the 183 GHz water transition. 



1. Introduction 

Water is one of the most important molecules in interstellar clouds in general and in star- 
forming regions in particular. In warm regions (T >100 K) close to the p rotostar, water is 



prom i nent with gas abundances up to 3x10 4 w.r.t. H ? , even higher than CO ( Cernicharo et al 



19901 Ivan Dishoeck & H elmich 



1996, Harwit et al 



1998). Besides the inner regions of 



protostellar envelopes these high abunda nces are also foun d where powerful jets from the 
protostar interact with the surroundings (INisini et al.l 1 19991) . In contrast, water abundances 
are as low as 10 -8 -10~ 9 in cold (T <100 K) envelope regions as evidenced by Infrared Space 
Observatory (ISO) a nd Submillimeter Wave Astronomy Satellite (SWAS) observations (e.g 
Boonman et al]|2003|) . Emission from water molecules is very difficult to observe, due to the 
limits imposed by the Earth's atmosphere. Isotopologues such as d euterated water can be 
observed from the ground (e.g. Schulz et al. 1991 . |Parise et al. 2005 ) The only transitions of 
the main water isotope that can be observed regularly are low-frequency maser transitions. 
Most famous is the water maser at 22.2 GHz, the 6i6-523 transition, regularly observed 
using radio telescopes and interferometers, such as the Very Large Array (VLA). In star 
forming regions with embedded sources of high luminosity (>100 h^) 22.2 GHz water maser 
emission is commonly dete cted and used to probe gas kinematics (e.g. lMoscadelli et al.l l2006. 
Goddi & Moscadellill2006|) . The emission is found to be variable on timescales of a day to 



Letter to the Editor 



2 



a month, probably related to variations in the accretion disk or the outflowing material (e.g. 
Pashchenko & Lekhtll2005l) . In regions with embedded sources of lower luminosity, 22.2 GHz 



water maser emission is dete cted less frequently t han in high-mass sources, but again with 



cons iderable variability (e.g. IWilking et al.l 1 19941 IClaussen et al.1 1 19961 : Il998l iFuruya et al, 



20031) . Imaging wit h the VLA detected the m aser emission within several hundred AU of low- 



mass protostars (e.g lFuruya et al.l 1 1 9991 ; 1200 If) . while Very Long Bas eline Array observation s 



indicate that the location may even be closer to the protostar (e.g iMoscade 



Both water and methanol masers have been observed to be related to disks (Torrelles et al 



li et al 



2006). 



1996, Moscadelli et al. 2006), while other observations of masers have been associated with 



outflows (e.g. IClaussen et al.l 1 1996L Furuya et al.l 120031) . However, the excitation conditions 
for the 22 GHz water maser l ine are very high density (>10 8 cm" 3 ) and temperature (2000 K> 
T >200 K) (|Yates et al.ll 1997b . a combination of conditions that is rare in low-mass protostars. 

Another water maser transition is found at 183.3 GHz. This 3i3-2 2 o transition has been 
used in extra-ga lactic and galactic studi es, primarily with the 30m telescope at Pico Veleta, 
Spain (e.g. ICernicharo et al.ll 19941 : 1 19960 . Statistical equilibrium calculations combined with 
the observation of extended emission conclude that relatively low temperatures (T~ 150 K) 
and densities (10^-10 6 c m" 3 ) can already invert the populations of the 183.3 GHz transition 
(ICernicharo et al] 1 19941) . In star forming regions it has been detected in several low- and 
high-mass protostellar sources such as Orion, W49N, H7-11 and L1448-mm. Most lines 
consist of a broad component, superposed with a strong narrow line, presumably a maser line. 
The broad component was found to be spatially extended thermal emission after having been 



mapped in some high-mass star form ing regions such as Orion and W49N (|Cernicharo et al. 



1 994 iGonzalez- Alfonso et al. 1995). Extended thermal emission was also found in low-mass 



sources such as HH7-1 1 (jCernicharo et al.ll 19961) . 

Serpens SMM1 (referred to as SMM1) is a low-mass Class sour ce in the Serpens 
cluster (D= 250 pc) and has been studied at (sub)millimeter wavelengths by lHogerheijde et aT 



(|1999j) . This source was selected as part of a pilot study to observe the 183 GHz line due to its 
inclusion in the Water in Star-forming regions with Herschel (WISH) program^ as a source 
where many water lines will be targetted with spatial resolutions of 9-40". It is relatively 
luminous ( L=20.7 L n ) an d thus an ideal target for water observations. It was observed with 
ISO-LWS (Larsson et alJl2002h . where numerous thermal water l ines were found in the large 
120" beam. It dri ves a powerful hi ghly collimated radio jet (|Curiel et al.l Il993|) . a large 
mol ecular outflow (IWhite et al.lll995|) and was included in the 22 GHz water maser survey 
by Furuya e t al. (20 031), wh e re sev eral maser spots were found. These masers were in turn 
studied by iMoscadelli et al] (I2006|) who found the masers originating within a circumstellar 
disk based on their position and proper motions. Observations with SWAS (beam = 3.3' by 
4.5') of the lio-loi water lin e found water assoc iated with the outflow with an abundance of 
o-H 2 of a few times 10" 7 (Franklin et al.ll2.Q08Q in both red and blue outflow lobes. 

In this letter, we present first interferometric observations using the SMABof the 183 

t See http://www.strw.leidenuniv.nl/WISH 

§ The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 



Letter to the Editor 



3 



GHz maser emission from a region of low mass star formation, associated with the protostar 
SMM1. 

2. Observations 

Observations were taken in September 2008 using the SMA in a compact configura- 
tion(maximum baselines ~75 m) under excellent weather conditions ( PWV < 1 mm). While 
the data set consisted of just one uninterrupted hour of usable data for the target, SMM1, 
(RA= 18:29:49.8, Dec=01: 15:20.5) it was enough to provide for unique spectral imaging of 
water masers. 

The SMA frequency coverage includes separation of the two sidebands, each covering 
2 GHz with an IF coverage of 4-6 GHz. The receiver was tuned to place the water line at 
183.310117 GHz in the lower sideband. At this frequency and configuration combination, 
the beamsize is 3.1" x 4.0". Due to pressure broadening, the opacity of atmospheric water 
vapor is spread in a roughly gaussian shape over several GHz, with a maximum at the rest 
frequency. To minimize the total average opacity across the full 2 GHz lower sideband (and 
thus increase sensitivity to continuum emission near the water line) the line was placed near 
one edge of the lower sideband, providing the best possible sensitivity at the other end of the 
sideband. 

The passband shape is a mix of instrumental bandpass characteristics and the atmospheric 
water absorption across the band. The instrumental terms of the atmospheric absorption were 
calibrated and removed by using high SNR observations of Jupiter. However, Jupiter was at a 
lower mean elevation than SMM1, and thus the atmospheric absorption was overcorrected for 
the target. Additional calibration of the atmospheric water vapor line shape was performed 
using spectrally smoothed data from 1751+096, a nearby bright quasar. Finally, under the 
assumption that an individual maser component is spatially unresolved, amplitude and phase 
self-calibration was performed using the brightest maser feature from SMM1 to improve the 
imaging reliability, though in the process we lose direct absolute position information. These 
absolute positions were retrieved by independently calibrating the upper sideband at 194 GHz 
using quasars and applying that solution to the already self-calibrated lower sideband. 

As the observations were done near the center of the atmospheric water line, we estimate 
the flux calibration error to be 40%, significantly larger than typical for observations well 
away from water lines. However, the errors in the relative fluxes are much smaller. 

For continuum purposes, additional observations were taken in July 2009 using the SMA 
in very extended mode with baselines of up to -400 kA. These observations were done at 220 
GHz to avoid the water absorption and data were taken for 5 hours. Data were calibrated on 
1751+096, with uranus and 3c454.4 as bandpass calibrator and Callisto as flux calibrator. 

the Academia Sinica. 
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3. Results 



3.1. Maser emission 

The emission of the 183 GHz line shows a main line that is very bright (460 Jy/beam) with 
very narrow line-width (1 km s _1 ) (Figure 1). Inspection of the map reveals two more line 
components, with brightnesses of 243 and 328 Jy/beam. If one assumes that the emission fills 
the beam, the brightness temperatures are 1871, 982 and 1325 K respectively. Brightness 
temperatures are even much higher if one assumes a smaller source. Typical brightness 
temperatures of 22 GHz maser emission can be on the order of 10,000 K from very small 
regions. However, without higher resolution observations, it is impossible to conclude if 
these 183 GHz spots originate in similarly small spots or from larger regions. Even assuming 
it fills the beam, such temperatures are much higher than the ambient kinetic temperatures 
expected in protostellar envelopes and are thus not thermally excited. 

The maser spots are labelled '1', '2' and '3' in Figured] Properties of the maser emission 
in each spot can be found in Table [Q while spectra at each spot can be found in Fig. [Q Several 
properties immediately jump out. First the center velocity of each of the spots decreases as a 
function of distance to the central protostar. Spot 3, which is the furthest away from the source, 
has the lowest velocity. Considering that the velocity of the Serpens SMM1 source itself is 
estimated at 8.7 km s" 1 from submillimete r single-dish and interfe rometric observations of 
molecular lines such as H 2 CO and N 2 H + (|Hogerheijde et al.lll999l) . this corresponds to the 
most blue-shifted velocity. Spot 1 is the only spot that is slightly red-shifted. However, the 
amount of both red- and blue-shifting is small. Second, the widths of the maser lines are not 
equal. The line associated with Spot 1 is a factor of 2 narrower than the other two, possibly 
indicating different excitation conditions. Third, the sp o ts see m to lie directly i n the path 
of the red- shifted outflow as determined by ICuriel et al.l (| 19931) and IWhite et al.l (| 19951) . In 
Curiel et al. Jl993|) it is also seen that at 3.6 cm continuum, small emission peaks appear at 
the red-shifted flow, but not the blue-shifted flow. These spots are at comparable positions as 
the maser emission (See also § 4) 



3.2. Thermal emission 

In the observations of Orion-KL, W49 and HH 7-1 1 dcernicharo et all 19941 : tl99~6l) . the 183.3 
GHz emission also showed a broad emission profile that was thermal in nature. Mapping 
indeed showed that such thermal emission existed over larger areas. In our interferometric 
observations, we do not detect any thermal emission, down to a level of 7 Jy equal to 16.2 
Kelvin, assuming it would fill the central beam. Due to the spatial filtering, it is possible that 
the emission is present at larger scales (>20"). However, this is not likely as in contrast with 
Orion, the interstellar radiation field ne ar SMM1 is not strong w ith no more massive stars 
present. If one adopts the models from Ivan Kempen et al] (I2008|) and simulates the results 
through the beam and spatial filtering of the SMA, one expects a line of 6 to 20 K within the 
central beam, equivalent to 2.5 to 8.4 Jy (depending on the combination of inner and outer 
water abundances, which range from 10 4 to 10 b for material warmer than 100 K and 10~ 6 
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to 10~ 8 for material colder than 100 K). The emission predicted by this model is sufficiently 
compact, spanning less than 5", that spatial filtering by the SMA would have no appreciable 
effect on a detection. 



3.3. Continuum 

In addition to observing the line, the continuum was extracted from the same observations 
using the upper sideband at 194 GHz, which was free of lines. The resulting map is shown 
in the left part of Fig. El with the amplitude versus UV distance plot shown below it. The 
peak flux of the map is 0.52 Jy/beam. Due to the limited UV-coverage, only 5 bins could be 
made. Overplotted is the expected visibility am plitude of a spherical model envelope based 
on SCUBA observations, modelled with Dusty dlvezic & Elitzur II 19971) and sampled with the 
194 GHz SMA coverage (Kristensen et al. in prep). This model has an envelope mass of 4.1 
solar masses. The emission at longer baselines clearly is not fitted by the envelope model. 
Crosses represent the residual visibilities. If one ignores the drop in emission at the shortest 
wavelengths, this represents a point source with a total emission of ~ 0.55 Jy. An upper limit 
to the disk radius of <325 A U can deduced from this fi t, which is in agreement with a compact 
source of ~90 AU found by iHogerheijde et al.1 (|1999). 



Subsequently, observations were done of the continuum in a very extended configuration 
(k/l=400) with an rms of ~2 mJy/beam. The resulting image is shown in the right side of the 
Fig. |2] The disk is resolved there, with a radius of ~100 AU. The plot of the amplitude versus 
UV distance shows that emission at baselines larger than 50 k/i remains, until about 250 k/l 
after which the sensit i vity b ecomes too low for imaging the inner regions of the resolved 
disk. Moscadelli et al.l (|2006) find the disk to be aligned in a NE-SW direction, but that is not 
confirmed . 

The integrated flux of the compact emission in a 6" radius is 2.3 Jy, with an 
associated mass of 2.7 M fi , assuming a temperature of 20 K and OH5 dust emissivities of 



Ossenkopf & Henningl (119941) . The total dust mass of 4.1 M Q also agrees well (Kristensen et 
al. in prep), as most of the cold dust e mission on scales l arger t han 20" is filtered out by the 
lack of baselines smaller than 10 k/l. IHogerheijde et al. d 1999b find a flux of 2.65 Jy at 1.4 
mm (214 GHz), which is consistent with a v 3 5 dependence. Note that their mass of 8.7 M is 
significantly higher due to the adopted distance of 400 pc. 

If we assume that all the emission in the very extended image, which equals to 0.21 Jy, is 
originating in a resolved disk, then the total mass of the disk i s ~0.1 M , less than l/2 th of the 
envelope mass. When compared to the sources observed in 1 Jprgensen et al.l (I2007|) . one can 
conclude that although the envelope around SMM1 may be exceptionally massive compared 
to the disk and the disk on the small side, the observed profiles and masses are not unusual for 
a Class source. 
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4. Discussion 



4.1. 183 GHz maser excitation 

Two main studies have been carried out where the excitation conditions of the 183 GHz 



maser have been investigated. Following lYates et al.l (I1997|) . the regime in which both 22 
GHz and 183 GHz transitions can mase is at similarly high densities of 10 8 to 10 10 cm -3 
and temperatures of 200 to 2000 K. Such spots would be constrained to very small regions, 
probably about 1 AU i n size at the maximum. However, a statistical equilibrium study of 
maser emission done by iCernicharo et al.1 (11994) found different excitation conditions for the 
183 GHz line, equivalent to densities of 10 5 -10 6 cm" 3 and temperatures in excess of 50 K. 
Such conditions are much more prevale nt near low-mass protostars, e.g. within the inner 
envelope regions, than those suggested by I Yates et al.l (119971) . For SMM1, the maser positions 
are a minimum of 675 AU away from the central protostar, a much larger distance than the 
size of the disk fou nd in the UV-amplitud e fitting (see above) and the location of the 22 
GHz emission from lMoscadelli et al.l (|2006|) . thus making it very unlikely that the two maser 
originate in the same region. The densities of the protostellar envelopes at the distances of the 
maser spots (+/- 700 to 1200 AU) are 10 6 cm" 3 ( J0rgensen et al.ll2002|) . However, envelope 
temperatures derived by low-excitation submillimeter lines are generally much lower (~50 
K), except in regions where the outflow and envelope interact. 



4.2. Maser origin 

From the emission profiles in Fig. [Q the inherent properties of the 183.3 GHz maser 
excitation, comparison with the 22 GHz maser as well as the range of properties listed in 
Tabled! the most likely origin for the emission of the 183 GHz masers is the cavity wall that 
the outflow jet creates by interacting with the spherical protostell ar envelope. This can indee d 
be an origin for the 22 GHz maser line in high-mass p rotostars dGoddi & Moscadelli 2006 ). 
Delamarter et al.l (|2000l) and ICunningham et al.l (|2005l) performed detailed hydrodynamical 
simulations that show that outflows can sculpt the protostellar envelope, creating high- 
density regions (>10 7 cm -3 ) in the cavity walls due to Kelvin-Helmholtz and Raleigh-Taylor 
instabilities. Temperatures at the cavity walls can reach a few hundred Kelvin from passing 
shoc ks and illumination by U V photons, an effect seen in recent observation of CO 6-5 and 
7-6 (|van Kempen et al.ll2009r ). Water will also be liberated from the grains, allowing for high 
abundances. This environment, created uniquely in the cavity walls, is able to provide the 
necessary conditions for the populations to invert, creating the maser spots. A major unknown 
influence is probably the path length of the emission through the medium. 

One problem however, is t he velocity field, a s the m aser sp ots seem to be infa l ling in stead 
of outflowing. As seen from Delamarter et al.l (|2000l) and ICunningham et al.l (|2005r ). the 
walls are very turbulent with infall and outflow material influencing each other in a Kelvin- 
Helmsholtz instability. It would pose little problem for some material in the cavity walls 
to be slightly infalling, explaining the presence of blue-shifted material in the red-shifted 
outflow. The changing path-lengths of the photons as they cavity walls evolve will also 
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produce variability in the spots. 
5. Conclusion 

The conclusions of this letter are as follow 

• The 183 GHz emission of Serpens SMM1 is masing from three distinct maser spots 
at 500-1200 AU, aligned in the direction of the red-shifted outflow, two of which are 
slightly blue-shifted. 

• The extended dust continuum emission detected with the SMA in Serpens SMM1 comes 
from a protostellar envelope with an estimated mass of 2.7 M . 

• A resolved component detected on baselines up to ~400 kA indicates the presence of disk 
with an estimated mass of 0.1 M . 

• It is theorized that the maser spots originate in Kelvin-Helmholtz instabilities in the 
red-shifted outflow, created by the jet interactions with the envelope. Such interactions 
would temporarily create small regions of high density and temperature, providing the 
conditions for the 183 GHz transition to mase in low-mass protostellar environment. 

• No compact or extended thermal water emission is detected, in agreement with a model 
of a passively heated envelope. 

Research into the 183 GHz maser is promising as an additional constraint to jet mechanics 
and interaction of the jet with its surrounding material, both the envelope and the outflow. In 
the next few years, the receivers of the SMA provide a unique opportunity to study the water 
maser emission at this wavelength in many low-mass YSOs. In the long run, ALMA will 
be able to observe these lines at higher sensitivity and resolution, possibly probing the scales 
at which these masers emit. Combined with high-resolution VLA and VLBI observations of 
the 22 GHz maser, such observations can truly probe the physical structure of maser-emitting 
regions. 
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Figure 1. Moment map of the three maser spots in SMM1. Maser spots are labeled with 
their respective number. An X at the center 0,0 position marks the maximum intensity of the 
continuum (see Fig. 01, which coincides with the location of the protostar. The spectra at the 
three maser spots are shown . Due to the proximity of the spot 2 to both spot 1 and 3 (~ 1 
beamsize), the spectrum there is a blend of all three maser spots. The beamsize is to 2.7" by 
4.0" with a P. A. of 49°. 



Table 1. Properties of the maser spots 



Maser spot 


FWHM 


^peak 




^LSR 


Spat. Offset 




km s _1 


Jy/beam 


K 


km s _1 


RA,Dec (",") 


1 


0.7 


463 


1871 


9.2 


-3.0,0.2 


2 


1.5 


243 


982 


7.4 


-4.5,2.2 


3 


1.3 


328 


1325 


6.3 


-6.7,3.8 



a assuming the emission fills the beam 
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Figure 2. TOP: The continuum map of SMM1 at 187 GHz. Levels are in 2cr, Act, 6cr,... with 
<x=15 mJy/beam. BOTTOM: The UV distance versus amplitude plot, with the mean visibility 
amplitude in that annulus on the vertical axis and baseline length on the horizontal. The 
triangles show the observed visibilities of SMM1 at 187 GHz. Errors shown are the standard 
deviation in the mean, while the solid line shows the expectation value assuming no signal. 
The dashed line shows the expected visibilities for an envelope model at 194 GHz (Kristensen 
et al. in prep). The crosses represent the residual visibilities if the envelope model is subtracted 
from the observed visibilities. 



